function[r]=f(n,k)
r=prod(1:n)/(prod(1:k)*prod(1:(n-k)));
end

